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ABSTRACT 

Turbulence with a large magnetic Reyonolds number, generically leads to rapidly 
growing magnetic noise over and above any mean field. We revisit the dynamics 
of this fluctuating field, in homogeneous, isotropic, helical turbulence. Assuming 
the turbulence to be Markovian, we first rederive, in a fairly transparent manner, 
the equation for the mean field, and corrected Fokker-Plank type equations for the 
magnetic correlations. In these equations, we also incorporate the effects of ambipolar 
drift which would obtain if the turbulent medium has a significant neutral component. 
We apply these equations to discuss a number of astrophysically interesting problems: 
(a) the small scale dynamo in galactic turbulence with a model Kolmogorov spectrum, 
incorporating the effect of ambipolar drift; (b) current helicity dynamics and the 
quasilinear corrections to the alpha effect; (c) growth of the current helicity and 
large-scale magnetic fields due to nonlinear effects. 



Subject headings: Magnetic fields; turbulence; dynamo processes; ambipolar drift; 
helicity; Galaxies 



1. Introduction 

The origin of large-scale cosmic magnetic fields remains at present, a challenging problem. In 
a standard paradigm, one invokes the dynamo action involving helical turbulence and rotational 
shear, to generate magnetic fields ordered on scales much larger than the turbulence scale (cf. 
Moffat 1978, Parker 1979, Zeldovich et al. 1983). However, turbulent motions, with a large enough 
magnetic Reynolds number (MRN henceforth), can also excite a small-scale dynamo, which 
exponentiates fields correlated on the tubulent eddy scale, at a rate much faster than the mean 
field growth rate (Kazantsev 1968, Zeldovich et al. 1983 and see below). A possible worry is that, 
these small-scale fields can come to equipartition with the turbulence much before the large-scale 
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field has grown appreciably and may then interfere with the large-scale dynamo action (cf. 
Kulsnid and Anderson, 1992 (KA)). Indeed the efficiency of the dynamo to produce the observed 
large-scale field has come under increasing scrutiny (Cattaneo and Vainshtein 1991, Vainshtein 
and Rosner 1991, Kulsrud and Anderson 1992, Brandenburg 1994, Field 1996, Blackman 1996, 
Chandran 1997, Subramanian 1997). It appears therefore that understanding large-scale dynamo 
action due to helical turbulence is also closely linked with understanding the dynamics of the 
small-scale fluctuating component of the magnetic field. This forms the prime motivation for the 
present work, where we revist the dynamics of the fluctuating field ab initio. 

Fluctuating field dynamics is best studied in terms of the dynamics of magnetic correlation 
functions (see below). Kazantsev (1968) derived the equations for the longitudinal correlations 
in homogeneous, isotropic , Markovian turbulence, without mean helicity, using a diagrammatic 
approach. Vainshtein and Kitchatinov (1986) incorporated the effects of helicity, and derived 
equations for both helical and longitudinal correlations. They did not however give details of the 
algebra. We therefore first rederive these equations here, in a fairly transparent manner. Our work 
generalises the Kazantsev equations to helical turbulence, and corrects a sign error of Vainshtein 
and Kitchatinov (1986) in the helical terms, which could be important for understanding the back 
reaction of the small-scale field on the alpha effect (see below). 

The above mentioned works on small-scale field dynamics were also purely kinematic. As the 
field grows one expects the back reaction due to the Lorentz force to become important. The full 
MHD problem, involving both the Euler and induction equation, presents a formidable challenge, 
which we do not take up here. However in order to get a first look at the effect of nonlinear ities, 
we consider a simpler nonlinear modification to the velocity field due to the generated magnetic 
field, as would obtain for example in a partially ionised medium. Note that the gas in galaxies, 
proto stellar disks and in the universe just after recombination are all likely to be only partially 
ionised. In such a medium, the Lorentz force on the charged component will cause a slippage 
between it and the neutrals. Its magnitude is determined by the balance between the Lorentz 
force and the ion-neutral collisions. This drift, called ambipolar drift (Mestel and Spitzer 1956, 
Draine 1986, Zweibel 1988), can be incorporated as a field dependent addition to the fluid velocity 
in the induction equation. This makes the induction equation nonlinear and provides us with a 
model nonlinear problem to study. Such a modification of the velocity field has also been used by 
Pouquet et al. (1976) and Zeldovich et al. (1983) (pg. 183) to discuss nonlinear modifications to 
the alpha effect. 

We incorporate the effects of this nonlinearity while deriving the evolution equations for both 
the mean field and fluctuating field correlations, in section 3 and 4. The presence of a nonlinear 
term due to ambipolar type drift in the induction equation implies that lower moments couple 
to higher order moments. Some form of closure has to be assumed. We adopt here a Gaussian 
closure. The mean field and the magnetic field correlations then obey a set of nonlinear partial 
differential equations with the nonlinearity appearing as time-dependent co-efficients involving 
the average properties of the fluctuating field itself. In the sections which follow we apply these 
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equations to study a number of astrophysically interesting problems. 

We begin in section 5 by examining small-scale dynamo action in the galactic context. This 
problem was studied by KA who looked at the evolution of the magnetic energy spectrum in 
wavenumber space, taking the limit of large k (small scales). In contrast to their treatment, the 
co-ordinate space approach adopted here, allows us to implement the boundary conditions for the 
magnetic correlation function, at both large and small-scales. Using the WKBJ approximation, 
we derive conditions for growth of small-scale fields in a model Kolmogorov turbulence and study 
the eigen functions for the longitudinal magnetic correlation (see below). Most earlier work also 
discussed small scale dynamo action for the case when the turbulent velocity has a single scale (cf. 
Kleeorin et al. 1986, Ruzmaikin et al. 1989 and references therein). Our work extends this to the 
context where multiple scales are present. 

The nonlinear effects of ambipolar drift, on the small-scale dynamo are considered in section 
6. The results of section 5 and 6, suggests a useful visualisation of the spatial distribution of the 
small-scale dynamo generated field (cf. Zeldovich et. al, 1983). We argue that in the galactic 
context, the magnetic field generated by small scale dynamo action, concentrates into thin 
(perhaps ropy) structures. One of the aims of our work here is to lay a framework for a companion 
paper (Subramanian, 1997 : Paper I). In Paper I we have built upon the results obtained in 
sections 5 and 6, to discuss in detail how the small-scale dynamo may saturate in the galactic 
context, in a manner which preserves large-scale dynamo action. 

Section 7 concentrates on the kinematic evolution of the helical component of the magnetic 
correlations. The study of helical magnetic correlations has been mostly neglected in the literature, 
since they do not drastically affect the small-scale dynamo. However the evolution of current 
helicity is important in determining the back reaction of small-scale fields on the a-effect in the 
standard dynamo equation (cf. Pouquet et al. 1976, Zeldovich et al. 1983, Gruzinov and Diamond 
1994, Bhattacharjee and Yuan 1995). Numerical simulations by Tao et al, 1993 indicated that 
the alpha effect is drastically decreased by the growing small-scale field (see however Brandenburg 
1994). We will examine this issue in terms of our approach. 

In section 8, we study the nonlinear effects due to the current helicity. We point out that 
the current helicity may increase temporarily due to its non linear coupling with the logitudinal 
magnetic correlation function. Wc also comment on the possibility of a non linear dynamo driven 
by the magnetic alpha effect and on the ordering of small scale fields due to relaxation and 
selective decay. The final section presents a summary of the results. 

2. Mathematical preliminaries 

In a partially ionised medium the magnetic field evolution is governed by the induction 
equation 

(dB/dt) = V X (vj X B - ?7V X B), (1) 
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where B is the magnetic field, Vj the velocity of the ionic component of the fluid and rj the ohmic 
resistivity. The ions experience the Lorentz force due to the magnetic field. This will cause them 
to drift with respect to the neutral component of the fluid. If the ion-neutral collisions are rapid 
enough, one can assume that the Lorentz force on the ions is balanced by their friction with the 
neutrals. Under this approximation the Euler equation for the ions reduces to : 

Pil^ini^i - V„) = PiUin^D = [(V X B) X B]/(47r), (2) 

where pi is the mass density of ions, the ion-neutral collision frequency and v„ the velocity of 
the neutral particles. We have also defined here v/5, the ambipolar drift velocity. 

Using the Euler equation for the ions and substituting for Vj, the induction equation becomes 
the nonlinear equation 
(9B 

— = V X [v„ X B + a[((V X B) X B) X B] - r?V X B] , (3) 
where we have defined 

^T^Pil'in 

We have derived the above equation as one which describes the effect of ambipolar drift. 
However, one can also view Eq. ^ as describing a model nonlinear problem, where the nonlinear 
effects of the Lorentz force are taken into account as simple modification of the velocity field. 
Such a phenemenological modification of the velocity field has in fact been used by Pouquet et al. 
(1976) and Zeldovich et al. (1983) (pg. 183) to discuss nonlinear modifications to the alpha effect. 

The velocity field v„ is taken to be prescribed independent of the magnetic field. We will 
assume v„ has a turbulent stochastic component over and above a smooth component vq, that 
is v„, = Vq + v^. Since is stochastic, Eq. (^) becomes a stochastic partial differential equation. 
Its solution depends on the statistical properties of the velocity field v^. 

We assume v^^ to be an isotropic, homogeneous, Gaussian random velocity field with zero 
mean. For simplicity, in this work, we also assume to have a delta function correlation in time 
(Markovian approximation) and its two point correlation to be specified as 

<v^{-^,t)v^{y,s)>=r\r)5{t-s) (5) 

with 
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r^{r) = T^r^[6'^ - (^)] + TllC-^) + Ce,,frf. (6) 

Here <> denotes averaging over an ensemble of the stochastic velocity field vt, r = |x — y|, 

= — and we have written T^^ (r) in the form appropriate for a statistically isotropic and 
homogeneous tensor (cf. Landau & Lifshitz 1987). T/^L(r) and T]\f]\f{r) are the longitudinal and 
transverse correlation functions for the velocity field while C(r) represents the helical part of the 
velocity correllations. If vt is assumed to be divergence free (which we do here), then 
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The stochastic Eq. (||) can now be converted into equations for the various moments of the 
magnetic field. To derive these equations we proceed as fohows: (see also Zeldovich et al. 1983, 
Chapter 8) 



3. Mean field evolution 

Let the magnetic field at an initial time say t = be B(x, 0). Then, at an infinitesimal time 
5t later, the field is given iteratively by 

B(x,(5t) = B(x,0) + (5t 7/V^B(x,0) + / dtV x Ea(x, t) + / dt dsV x Eb(x, t, s) (8) 

JO Jo Jo 

where we have defined Ea(x, t) = V(x, t) xB(x, 0), Eb(x, t, s) = vt(x, t) x [Vx (vr(x, s) xB(x, 0))] 
and V = v„ + V£). We have also retained only terms which are potentially of order 5t, and so will 
survive in the limit 5t — > 0. Note that the last term in the above equation is potentially first order 
in 6t, because of the presence of the stocahstic turbulent velocity field v^. 

The time evolution of the mean magnetic field Bq =< B > can be deduced by taking the 
ensemble average of Eq. (|^) . As before <> denotes averaging over an ensemble of the stochastic 
velocity field v^-. On using the fact that vj- at time t is not correlated with either the initial 
magnetic field B(x, 0) or the initial perturbed field 5B(x, 0) = B— < B > and taking the limit 
6t ^ we get after some straight forward algebra 

^ = V X [vo X Bo + 2C7(0)Bo - (r/ + Tll{0))V x Bo]+ < V x v^j x B > (9) 
ot 

The effect of the turbulent velocity is to introduce the standard extra terms representing the 
a-effect with 

a = 2C(0) = --j< vt(x, t).V X vt(x, s) > ds (10) 
and an extra turbulent contribution to the diffusion 

VT = Tll(O) = IJ < vt(x, t).VT(x, s) > ds. (11) 

Over and above these terms the effect of ambipolar drift is to introduce an extra EMF 
represented by the last term in Eq. ^ . This extra term involves the third moment of the 
magnetic field. Similarly as we will see below a consideration of the equation for the magnetic 
correlation function will introduce the fourth moment of the magnetic field. This will lead to an 
infinite hierarchy of equations for the moments which can be only truncated by assuming some 
form of closure. We assume below, for analytic tractability, that the magnetic field fiuctuation 
5B = B — Bq is also a homogeneous, isotropic, Gaussian random field with zero mean. Its equal 
time two point correlation is given by 

< SB'{^,t)5B^y,t) >= M'\r,t) = Mn[6'' - (^)] + Ml(-^) + Hei^fr^ . (12) 
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(Here the averaging is a double ensemble average over both the stochastic velocity and stochastic 
5B fields, although we indicate only one angular bracket). The functions ML{r, t) and M]y{r, t) are 
the longitudinal and transverse correlation functions for the magnetic field while H{r, t) represents 
the helical part of the correlations. Note that H is proportional to the current-field correlation, 
rather than the field-vector potential correlation (see below). Since V.B = 0, Mn and Ml are 
related by 

The nonlinear term due to ambipolar drift, <Vxv£)XB>in Eq. (^) is then given by 

V X < 6B.V x6B> +a(Bo.V x Bo)]Bo - < dB.6B > +a(Bg)]V x Bq. (14) 
Using the form for the magnetic correlation function we have 

< 6B.V X 6B >= -6/7(0, t) and < 6B.6B >= 3ML(0,i). (15) 
So the mean magnetic field satisfies the equation 

V X [vo X Bo + ae//Bo -??e//V X Bo]. (16) 



dBo 



dt 
where 

aeff = 2C(0) - 4aH{0, t) + a(Bo.V x Bq) (17) 
Veff = r] + Tll{0) + 2aML{0,t) + aBg (18) 

The effect of ambipolar drift (or the field dependent addition to the fluid velocity) on the 
evolution of the mean field is therefore to modify the o-effect and the diffusion of the mean field. 
When one starts from small seed fields, the additional nonlinear terms depending on the mean 
field itself are subdominant to the terms depending on the fluctuating field for most part of the 
evolution; since we will find below that the small scale fields grow much more rapidly compared to 
the large-scale mean field. When ambipolar drift is taken into account, the small-scale fluctuating 
fields contribute an extra diffusion term to the mean field evolution, proportional to their energy 
density. Also the alpha effect is modified by the addition of a term proportional to the mean 
current aligned component (or current helicity, H(0,t)) of the magnetic field fluctuations. Some 
aspects of mean field dynamos incorporating ambipolar drift has been discussed by Zweibel (1988) 
and Proctor and Zweibel (1992). 



4. Evolution of the correlation tensor of magnetic fluctuations 



The derivation of these equations involves straightforward but rather tedious algebra. We 
therefore only outline the steps and the approximations below leaving out most of the algebraic 
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details. We start by noting that 

{dMij/dt) = {d/dt){<dB,{x,t)6Bj{y,t)>) 

= [{d/dt){<BiBj>)-{d/dt){<Bi><Bj>)]. (19) 

The second term in the square brackets is easy to evaluate using the equation for the mean field. 
The first term can be evaluated using Eq. @ and the fact that 

{d/dt){Bi{^,t)B,{y,t)) = Bi{^,t){dBjiy,t)/dt) + {dB^{^,t)/dt)Bj{y,t). (20) 

The resulting equation can again be solved iteratively to get an equation for (dMij/dt) which 
depends on the the turbulent velocity correlations Tij, the mean velocity field vg and the mean 
magnetic field Bq and most importantly a non-liner term incorporating the effects of ambipolar 
drift. We get 



dMjj 
dt 



< / 



+ < / 



+ < 



ds > 
ds > 



v^y,t) '^Ram{v'T{^,s)[Mrnq + B^{^)B^Q{y)) 
,p, ^v^{^,t) yRjim{vUy,s)[Mgm + 5g(x)S^(y)) 

yRjpg (<(y,t) yRglm{vT{y,s)Mim)) ds > 
~l~ ^ / Ripq 

(v^x,t) ''Rgin,{virix,s)Mmj)) ds > 

+ 7j[VlMi, + VlMij] +y Rjpg {vP{y)M,g) +^ R,pg {vP{i^)Mgj) 

+ yRjpg (< vl{y)6Bi{^)Bg{y) >) +^ R^pg (< T;^(x)5,(x)<5i?, (y) >) (21) 

where we have defined the operators 

""Ripg = eumempqid/dx'-) and ^Ripg = tiim^mpq{d / dy'') . (22) 



The first two terms on the RHS of Eq. ( pi]) represent the effect of velocity correlations on 
the magnetic fluctuations (Mjj) and the mean field {Bq). The next two terms the "turbulent 
transport" of the magnetic fluctuations by the turbulent velocity, the 5th and 6th terms the 
" microscopic diffusion" . The 7th and 8th terms the transport of the magnetic fluctuations by the 
mean velocity. The last two nonlinear terms give the effects of the back reaction due to ambipolar 
drift on the magnetic fluctuations. 

We note that the effects of the mean velocity and magnetic fields are generally subdominant 
to the effects of the fluctuating velocity and magnetic flelds. First the mean flelds vary in general 
on a scale much larger than the fluctuating flelds. If one neglects the effect of velocity shear due 
to the mean velocity, on the fluctating fleld, then one can transform away the mean velocity by 
going to a different reference frame. The mean magnetic fleld, ofcourse cannot be transformed 
away. But we will see that, it grows at a rate much slower than the fluctuating fleld. So when one 
starts from small seed magnetic fields, for most part of the evolution, we can neglect its effects. In 
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what follows we will drop the term involving the mean fields. Due to the above reasons, we also 
continue to treat the statistical properties of the magnetic fluctuations as being homogeneous and 
isotropic, and use Mjj(x,y, = Mij{r,t) as before. 

All the terms in the above equation, can be further simplified by using the properties of 
the magnetic and velocity correlation functions. In order to obtain equations for Ml and H, we 
multiply Eq. (^) by [r^r^^lr"^ and f-ijfrf and use the identities 

ML{r) = M,j{r\'^ It"), H{r) = Mije^jfr^ /{2r^). (23) 

We have given some steps in simplifying the first two terms in appendix A. The 3rd and 4th terms 
add to give a contribution 4:C{0)ejgm{dMim/dr'') + ITLLiO)"^"^ Mij to the RHS of Eq. (||) , hence 
justifying their being called "turbulent transport" of Mij. 

The last two nonlinear terms give the effects of the back reaction due to ambipolar drift on 
the magnetic fluctuations. In evaluating this term, we will neglect the effects of the subdominant 
mean field compared to the back reaction effects of the fluctuating field. In this case the nonlinear 
terms add to give a contribution -8aH{0,t)ejqjn{dMim/dr'i) + 4aML{0,t)V^Mij to the RHS of 
Eq. (|2^) . The Gaussian assumption of the magnetic correlations results in the nonlinearity of 
this term appearing as a nonlinearity in the coefficient, rather than the correlation function itself. 
Gathering together all the terms, we get for the coupled evolution equations for Ml and H : 



where we have defined 



KN = r] + TLLiO)- Tll {r) + 2aMi (0, i ) 
UN = 2C(0) - 2C(r) - 4ai/(0, t) 

(26) 



G 



dr r dr 



These equations together with Eq. (IC) for the mean magnetic field are an important result 
of this work. They form a closed set of nonlinear partial differential equations for the evolution 
of both the mean magnetic field and the magnetic fluctuations, incorporating the back reaction 
effects of ambipolar drift (or a magnetic field dependent addition to the velocity). For non-helical 
turbulence the equation for Ml excluding nonlinear effects was first derived by Kazantsev (1968). 
We note that Eq. (^) and Eq. (^) without the inclusion of the non linear terms due to ambipolar 
drift, have been derived in a different fashion by Vainshtein and Kichatinov (1986). We get exactly 
their equation (27) for Ml and H, except for a sign difference in front of the terms. We believe 
that our equations have the correct sign (see also below). 



- 9 - 



The terms involving kjv in equations ( p4D and (|25| ) represent the effects of diffusion on the 
magnetic correlations. The diffusion coefficient includes the effects of microscopic diffusion (rj) and 
a scale-dependent turbulent diffusion (Tll(O) — Tii{r)). The effect of ambipolar drift, under our 
approximation of Gaussian magnetic correlations, is to add to the diffusion coefficient an amount 
2aML(0,t) ; a term proportional to the energy density in the fluctuating fields. Similarly un 
represents first a scale dependent a-effect (2C(0) — 2C(r)) and the effect of ambipolar drift is to 
decrease this by 4ai?(0, t), an amount proportional to the mean current helicity of the magnetic 
fluctuations. Ambipolar drift has therefore very similar effect on the magnetic fluctuations as on 
the mean field. The addition of these terms makes equations (^) and (^) nonlinear, with the 
non linearity appearing in the coefficients. The term proportional to G{r), allows for for the rapid 
generation of magnetic fluctuations by velocity shear and the existence of a small-scale dynamo 
independent of the large-scale field (cf. also KA, Vainshtein and Kichatinov 1986). 

An important facet of MHD equations is the conservation of magnetic helicity, in the absense 
of microscopic diffusion. Magnetic helicity is defined as 

Im = J A.B d^x (27) 

where A is the vector potential. One can show, (cf. Moffat 1978) Im is conserved in the ideal 
limit. In case we have no mean fields, we also have Im = V < A.B > and 

1^ = -2r? < B.(V X B) > . (28) 

(Note that we are assuming that integration over a large volume V is same as the ensemble 
average). It is interesting to display this conservation explicitely from our equations. This will 
also act as check on some part of the algebra. 

For this let us go over to the vector potential representation of the fluctuating field and write 
5B = V X SA. We also define the equal time correlation fuction of the fluctuating component of 
the vector potential as P*-' =< 5A^5A^ >. Since we have taken the fluctuting magnetic field as 
a homogeneous, isotropic, Gaussian random field, Pij will also satisfy this property. So one can 
write in general 

P'\r,t) = Pn[5'' - {—)]+PL{—)+PHe^Jfrf. (29) 

Here the functions PL{r,t) and Pj\f{r,t) are the longitudinal and transverse correlation functions 
for the vector potential. PH{r,t) represents the helical part of the correlations and is related to 
the magnetic helicity of the fluctuating field, with 

PH{0,t) = - < 6A.5B > /6. (30) 



So to determine the evolution of the average magnetic helicity of the fluctuating field, we have to 
evaluate the time evolution of Pff(0, t). 
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One can easily relate correlation function H, representing the current-field correlation, to Ph, 
by using the definitions of the two quantities. We have 



Hir,t) 



2^.2 ^ij f^ilm^jpq^rnq,lp 



dr 



dr J 



Substituting this relation in the evolution equation ( |25| ) for H, we get 



0. 



Since all the correlations die off at spatial infinity the integral of this equation then gives, 

dPH(.r,t) 



dt 



-{2knH + aNML) 



(31) 



(32) 



(33) 



Now suppose we look at the evolution of the mean helicity of the fluctuating fields, in the 
absense of the mean magnetic field. Then conservation of Im implies conservation of PH{0,t). To 
see if this conservation obtains, take the r ^ limit of the Eq. (|3^). We get 

'-^ = -2,Hm (34) 

So when r/ = 0, the above equation shows that Pf/(0, t) = constant, independent of time, as 
required. Also when rj ^ 0, and the mean fields are zero, then Eq. (^4|) is identical to Eq. ( |28| ) , as 
required. It is heartening to note that the nonlinear equation that we have derived for the helical 
part of the magnetic correlations, incorporating the effects of ambipolar type drifts, does indeed 
embody magnetic helicity conservation in the appropriate limit, as required. 

Note that if the signs in front of the term had been different from the one we get, then 
the ambipolar drift terms in the equation for Ph would not have cancelled out when we take the 
r — > limit in Eq. (^). And we could not recover magnetic helicity conservation as above. This 
provides another check that we have indeed got the relative sign of these terms right. Also from 
the same argument, we note that one cannot have non-zero additions to the alpha term in the 
r — > limit (the — 4a-ff(0,t) term in a^) without a corresponding addition to the diffusion term 
(the 2aML{0,t) term in kn)- 

We now consider me of the implications of these equations for the evolution of small-scale 
fields. 



5. Kinematic evolution of A^l and the small-scale dynamo 

We begin by first studying the kinematic evolution of Ml, ignoring the effects of the nonlinear 
coupling terms due to ambipolar drift. In particular we consider the small-scale dynamo in galactic 
turbulence with a model Kolmogorov spectrum. The results obtained in this section will be used 
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to set the framework for Paper I. It also extends the analysis of the small-scale dynamo, due to 
turbulence with a single scale (cf. Zeldovich et al. (1983), Kleeorin et al. (1986)) to the situation 
when a range of scales are present. 

First we note that in many contexts the coupling term a^H, due to helicity fluctuations 
of the velocity, has negligible influence on the evolution of M^. A canonical estimate in the 
galactic context (see Zeldovich et al. 1983) is un ~ 2C(0) ~ {V'^t /h){VLT) ~ QrViL/h) and 
kn ~ Tll{Q) ~ V'^T ~ VL. Here where h is the disk scale height, O the rotation frequency, V , 
T and L are the velocity, correlation time, and correlation lengths for the energy carrying eddies 
of the turbulence. Assuming H ~ Mi/h, the importance of the coupling term a^H in Eq. (p^) 
to the other terms is ~ {un H / {i<-n Ml / L?)) ~ ^}T{L/h)'^ << 1 in general. (Even for r ~ L/V 
we generally have Or << 1). So it is an excellent approximation to neglect the coupling to H in 
examining evolution of Ml- We will refine this estimate in the next section after studying the 
evolution of H and see that the above approximation is even better. 

The evolution for Ml can then be transformed into a Schrodinger-type equation by defining 
$ = P^/k^Ml- We get 

where for a divergence free velocity field, the "potential" 

Uir,t)=n, + lT^^, + f-^-!^ + '-^. (36) 

The boundary condition on ML{r,t) is that it be regular at the origin and that ML{r,t) ^ as 
r — > oo. 

In the kinematic limit, note that k^V) Tll and hence the potential U are time independent. 



Equation (35) then admits eigenmode solutions of the form ^{r,t) = exp(2rt)<I>(r) where 

KN{d^^/dr^) - (T + U)^ = 0. (37) 

So there exists a possibility of growing modes with F > 0, if one can have U sufficiently negative 
in some range of r. The problem of having a small-scale dynamo and rapid growth of magnetic 
fluctuations, reduces to that of having bound states in the potential U. 

To see what this requires let us consider a model problem where the behaviour TLL{f) 
simulates Kolmogorov turbulence (Vainshtein 1982); 

TLLir) = ^[1-RV\j)^] forO<r<Z, 

= for/e<r<L 

= for r > L (38) 

Here Ic ~ LRe is the cut off scale of the turbulence, where Re = VL/v is the fluid Reynolds 
number and u is the kinematic viscosity. For Kolmogorov turbulence, the eddy velocity at any 
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scale I, is vi (X l^^^, in the inertial range. So the scale dependent diffusion coefficient scales as 
vil oc l^^^. This scaling, also referred to as Richardsan's law, is the motivation for the form of the 
scaling of T^^l with r which we have adopted for the inertial range of the turbulence. Note that 
the structure function Tll{i") must satisfy the condition T'j^x^{0) = 0, at the origin. So for scales 
smaller than Ic we have continued Tll from its value at r = Ic to zero, satisfying this constraint. 
(One can adopt smoother continuations for TiL{r) at r = Ic and r = L, in order to make the 
potential U continuous at these points. But this has little effect on the conclusions below, since $ 
is determined by integrals over U). 

The potential is then U = {2r]/r'^) as r ^ and U = 2{rj + riT)/r'^ for r > L. For Ic < r < L 
we have 

"_8(I)-2/3 _ (4/9)(r/L)V3 ^ 



(39) 



ri' (3/i2^ + (r/L)4/3) 

where Rm = (VL/rj) is the magnetic reynolds number at the outer scale of the turbulence. The 
Spitzer value of the resistivity gives r] = cm s . For numerical estimates we 

generally take V = lOkms^^ and L = lOOpc. For these turbulence parameters Rm = 3 x 10^^. 

The value of the potential at r = L is U ^ {V/L)[{2/Rm) - (4/9)]. If Rm < 9/2 one can easily 
see from the above expressions that U remains positive for all r. Note that for a bound state to 
obtain and for the small-scale fields to grow, one must have U negative for some range of r. So a 
neccessary condition for the small scale fields to grow is Rm >> 1. In fact the exact critical MRN, 
say Rc, for growth has to be determined numerically but one typically gets Rc ~ 60 (cf., Zeldovich 
et al 1983, Novikov et al. 1983; see below for a WKBJ derivation of this limit) 

Since the velocity at any scale /, say vi oc l^^^ for Kolmogorov turbulence, the MRN associated 
with eddies of scale is Rm{l) = vil/f] = Rm{l / L)^/^ . Using this one can also rewrite the potential 
U as 



^ 3/ 



8 ,r _ (4/9)(r/0^/3 6 L\ 

ri> {3/Rmil) + (r//)4/3) ^ R^{1) ^ r2 ■ 



(40) 



This is exactly of the SciiiiG form BS Eq. except tlicit ctnci Rrn 

are replaced hy l,vi and Rmi}) 
respectively. So a number of conclusions about the generation of small-scale fields can be scaled 
to apply to an arbitrary scale, /, provided we use the corresponding velocity scale vi and Reynolds 
number Rm{l) appropriate to the scale /. For example, the condition for excitation of small-scale 
dynamo modes which are concentrated at a scale /, is also Rmi}) = Rc » 1- Of course, exact 
scalability of the results only obtains well in the inertial range, far away from Ic and L. The MRN 
associated with eddies at the cut off scale is Rm{lc) = Vclc/v = V{lc/L)^^^L{lc/L)/ri = Rm/Re- 
(Here Vc is the eddy velocity at the cut-off scale). So if Rm/Re » Ij a potential well with U 
negative, extends up to the cut off scale of the turbulence, and these modes will also grow. 

We will be considering a largely neutral galactic gas and for this v is dominated by the neutral 
contribution. We take the neutral- neutral collision to be dominated by H-H collisons with a cross 
section an-H ~ 10~^^cm~^, leading to a kinematic viscosity u ~ vthi^/nHf^H-H)- For a thermal 
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velocity vth ~ Wkms and a neutral hydrogen number density nn ~ 1cm ^, as say appropriate 
for a young galaxy, we have u ~ lO^^cm^s"^, and 

i?, = — « 3 X lO^ioLioo (41) 

V 

where Vio = (y/10kms~^) and Lioo = (-L/lOOpc). So the condition Rm/Re > Rc in fact holds in 
the galactic context. 

It should also be noted that the value of the potential at any /, U{1) ~ vi/l, is the inverse of 
the turnover time of the eddies of scale /. The depth of the potential well at some scale I reflects 
the growth rate of modes concentrated around that I. And the growth rate of a mode extending up 
to r ~ say Ti ~ U{1) ~ vi/l ~ decreases with increasing I. So when Rm{lc) = Rm/Re > Rc, 

the small-scale fields tangled at the cut off scale grow more rapidly than any of the larger scale 
modes. These results have also been found in a different fashion by KA. To illustrate some of 
these points in a more quantitative fashion we have given a detailed a WKBJ analysis of Eq. (|37| ) 
in appendix B. Below we summarise the main results of this analysis: 

• The WKBJ analysis finds a critical value of the MRN, Rm = Rc ^ 60, for the excitation 
of the small-scale dynamo. Above this critical MRN the small-scale dynamo can lead to 
an exponential growth of the fluctuating field correlated on a scale L. We refer to the 
eigenmode which is excited for Rm = Rc as the marginal mode. Further, the equations 
determining Rc are the same if we replace {L,Rm) by {l,Rm{l))- Therefore, the critical 
MRN for excitation of a mode concentrated around r ~ / is also Rm{l) = Rc, as expected 
from the scale invariance in the inertial range. 

In the galactic context Rm » Rc', in fact, one also has Rm{lc) = Rm/Re » 1- Hence, 
small-scale dynamo action excites modes correlated on all scales from the cut-off scale Ic to the 
external scale L of the turbulence. 

• As expected, due to small-scale dynamo action, the fluctuating field, tangled on a scale 

grows exponentially on the corresponding eddy turnover time scale, with a growth rate 
r« ~ vi/l oc In the galactic context, with Rm{lc) = Rm/Re » Rc, the small-scale 

fields tangled at the cutoff scale grow more rapidly than any of the large-scale modes. 

• The WKBJ analysis gives a growth rate Fc = (?;c/Zc)[5/4 — co(ln(i?m/-Re))~^] with cq = 7r^/12 
for the fastest growing mode. Note that this is only weakly (logarithmically) dependent on 
Rm, provided Rm is large enough. 

To examine the spatial structure for various eigenmodes of the small-scale dynamo, it is more 
instructive to consider the function w{r,t) =< (5B(x, t).(5B(y, t) >, which measures the correlated 
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dot product of the fluctating field {w{0) =< (5B^ > ). Firstly there is a general constraint that can 
be placed on w{r). Since the fluctuating field is divergence free, we have 



since Ml is regular at the origin and vanishes faster than r~"^ as r ^ oo. Therefore the curve 
r'^w{r) should have zero area under it. Since w{0) =< (5B)^ >, w is positive near the origin. And 
the fluctuating field points in the same direction for small separation. As one goes to larger values 
of r, there must then values of r, say r ^ d, where w{r) becomes negative. For such values of r, 
the field at the origin and at a separation d are pointing in opposite directions on the average. 
This can be intepreted as indicating that the field lines, on the average are curved on the scale d. 

• We find that, in the case Rm/Re » 1) w{r) is strongly peaked within a region 

r = « lc{Rm/Re)^^^^ about the origin, for all the modes. Note that is the diffusive 
scale satisfying the condition ry/r^ ~ Vc/lc- For the most rapidly growing mode, w{r) changes 
sign across r Ic and rapidly decays with increasing r/Zc- For slower growing modes, with 
Ti ^ vi/l, w{r) extends up to r ~ Z after which it decays exponentially. 

• For the marginal mode w{r) peaks within a radius r ~ L/Rc^^, changes sign to become 
negative at r ~ L and dies rapidly for larger r/L. 

We should point out that a detailed analysis of the eigenfunctions can be found in Kleeorin 
et. al. (1986), for the simple case when the longitudinal velocity correlation function has only a 
single scale. Their analysis is also applicable to the mode near the cut-off scale of Kolmogorov 
type turbulence. These authors also give a pictorial intepretation of the correlation function, in 
terms of the Zeldovich rope-dynamo (cf. Zeldovich et al. 1983). If one adopts this intepretation, 
the small-scale field can be thought as being concentrated in rope like structures with thickness 
and curved on a scale upto ~ I for a mode extending to r L In Paper I, we also elaborate 
on a qualitative picture of the mechanism for the dynamo growth of small scale fields, and the 
generation of ropy fields from general initial conditions. 

As the small sclale fields grow the back reaction due to ambipolar drift will become important. 
We now turn to the nonlinear effects on the small scale dynamo to see when this can lead to 
small-scale dynamo saturation. 



w{r, t) 



1 r 3,. 
^— r^ML 
dr 



(42) 



so 




(43) 



6. 



Non-linear effects on the small-scale dynamo 



As the small-scale fields grow, ambipolar drift adds to the diffusion coefficient an amount 
2aMi(0, t) ; a term proportional to the energy density in the fluctuating fields. Similarly, 
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ambipolar drift leads to a decrease of by AaH{Q,t), an amount proportional to the mean 
current helicity of the magnetic fluctuations. When considering the evolution of the longitudinal 
correlation function, we once again neglect the subdominant effect of the the coupling term aj^H, 
due helicity fluctuations (sec above). 

As Mi(0,i) grows, its effect then is simply to change r] to an effective 

riambi=V + '^aML{0,t) (44) 

in the expression for the potential U{r,t). One can define an effective MRN, for fluid motion on 
any scale of the turbulence 

"^""^"^^^^ = ^ " 2aMU0,t) ^^^^ 

where vi = {l/Ly^^V as before. 

As the energy density in the fluctuating field increases Rambii^) decreases. Firstly, this makes 

1/2 

it easier for the field energy to reach the diffusive scales ~ ^c/ Rambi^^c)^ from a general initial 
configuration. More importantly as this happens the potential well disappears, first at small-scales 
and then progressively at larger and larger scales. This means that Ml will grow slower and 
slower. The detailed evolution of will be complicated. However, the decrease of Rambi suggests 
one possible nonlinear saturation mechanism for the small-scale field. The possibility that the 
system finds the stationary state with (dM^/dt) = 0. In such a state. Ml is independent of time. 
So, the condition on the critical MRN for the stationary state to be reached, will be identical to 
that obtained in the kinematic stage. Prom the discussion on the kinematic evolution of Ml, the 
stationary state is an eigenmodc for the system which obtains when the energy density of magnetic 
fluctuations has grown such that 

VL 

RambiiL) = T-r-, — r = Rc- (46) 

amOiK ) n + 2aML{0,t) ^ ' 

So if Rambi{Li) decreases to a value i?c ~ 60, dynamo action will stop completely. 

Let us consider now whether this condition can obtain in turbulent, partially ionised galactic 
gas. Take for example V ~ 10km s~^ and L ~ lOOpc appropriate for galactic turbulence. Also let 
us assume that the galaxy had very nearly primordial composition in its early stage of evolution: 
then the ions are mostly just protons and the neutrals are mostly hydrogen atoms. We estimate 
in Paper I that PiVin = riiPn < crv >eff with < av >e//~ 4 x 10~^cm^s~^. Here rii is the ion 
number density. We then have 

where f{l) = Bf /{ATrpnvf) is the ratio of the local magnetic energy density of a flux rope curved 
on scale I, to the turbulent energy density p^vf /2 associated with eddies of scale I. Using the 
value of as determined above and putting in numerical values we get 

Q{1) = ~ 1.8 X 10V_2(|)2/3LiooV-io' (48) 
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where n_2 = (rij/lO ^cm and we have assumed a Kolmogorov scaUng for the turbulent 
velocity fluctuations. 



turbulence, the MRN incorporating ambipolar drift is likely to remain much larger than Rc for 
most scales of the turbulence, even when the field energy density becomes comparable to the 
equipartition value. Only if the galactic gas is very weakly ionised, with nj < 10~^'^cm~^, can 
ambipolar drift by itself lead to small-scale dynamo saturation. Such small ion densities may 
perhaps obtain in the first collapsed objects in the universe, which collapse at relatively high 
redshifts (cf. Tegmark et al. 1997). However for most regions of a disk galaxy, the ion density is 
likely to be larger. So one expects the field to continue to grow rapidly, even taking into account 
ambipolar drift. Note also that the growth rates for the small-scale dynamo generally depends 
only weakly on the MRN, provided the MRN is much larger than R^. (see section 5, and Kleeorin 
et al. 1986). Therefore, we still expect the small-scale dynamo-generated field to grow almost 
exponentially on the eddy turn around time scale, as long as Rambi >> Rc- The spatial structure 
of the fluctuating field will also remain ropy, as in the kinematic regime, as long as Rambi >> be- 
llow then does the small-scale dynamo saturate in galaxies? In Paper I we consider this 
question in some detail by examining other nonlinear feedback processes which could limit the 
growth of the galactic small-scale dynamo, taking into account also ambipolar drift. We briefly 
suumarise our findings here for the sake of completeness. The reader is referred to Paper I for 
details. 

It turns out that the effect of the growing magnetic tension associated with the small-scale 
dynamo generated field is crucial. This tension acts to straighten out the curved flux ropes, 
while frictional drag damps the magnetic energy associated with the wrinkle in the rope. Also, 
small-scale flux loops can collapse and disappear. For a significantly neutral geas, these non-local 
effects turn out to operate on the eddy turnover time scale, when the peak field in a flux rope has 
grown to a few times the equipartition value. Their net effect is to make the random stretching 
needed for the small-scale dynamo inefficient and hence saturate the small-scale dynamo. However, 
the average energy density in the saturated small-scale field is sub equipartition, since it does not 
fill the volume. It is probable then that the small-scale dynamo generated fields do not drain 
significant energy from the turbulence, nor convert eddy motion of the turbulence on the outer 
scale to wave-like motion. So the diffusive effects needed for the large-scale dynamo operation can 
be preserved. This picture of small-scale dynamo saturation obtains only when the ion density is 
less than a critical value of n? ~ 0.06 — 0.5cm~^(n„/cm~^)^'''^. 

We note in passing that in the case of nearly neutral disks around protostars, the parameters 
could indeed be such that small-scale dynamo action saturates due to purely ambipolar drift (cf. 
Proctor and Zweibel 1992). We now turn to the the evolution of the current-field correlations. 



One can 




typical parameters associated with galactic 
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7. Kinematic evolution of current helicity and a-effect suppression 



The evolution of H is governed by Eq. ( pq ) derived in section 4. We consider first the 
kinematic linear evolution when the nonlinear terms in kn and aN are ignored. We saw above 
that the coupling of H to Ml is unimportant for the evolution of the longitudinal correlation 
function Ml, which could grow due purely to self-coupling terms. However for the evolution of H, 
its coupling to M^ provided by the ajy term cannot be neglected, because without the forcing by 



the Ml term in Eq. (25), any initial distribution of H decays with time due to diffusion. To see 
this multiply (25) by 2Hr'^KN and integrate over all r, neglecting the a^Mi coupling term. We 



get 



d_ 

di 



J H^KNr^dr =2KNHr^ — {2KNH)\^ - J drr^[— {2k NH)f . (49) 

Neglecting the boundary term for a sufficiently rapidly falling H we see that any smooth H will 
decay with time in the absence of the coupling to M^. 

So we have to consider the full inhomogeneous partial differential equation (^) for taking 
Ml to be a given function of r and t determined by the analysis of section 5. The solution to ( p5|) 
can be written then as the sum of the solution to the homogeneous equation Hq and a particular 
solution, say Hp, forced by the presence of Ml- Hq decays with time, as shown above and therefore 
the asymptotic evolution of H in time will be governed by Hp. To determine Hp{r, t) we use the 
Green's function technique (cf. Burton 1989). We are particularly interested in determining Hp 
when the longitudinal correlation function Ml = exp(r„i)M„(r), where the function Mn{r) is an 
eigenmode of section 5, with growth rate r„. In this case we can write Hp{r,t) = Hn{r) exp{rnt), 
with Hn{r) satisfying the time independent differential equation 

TnHn{r) = li^ [r^^{2KNHn + aNML)^ (50) 

The boundary conditions are that Hn{r) is regular at r = and Hn{r) ^ as r ^ cxd, as before. 



For general KAr(r) and a^ir) there are no general solutions to (|50|). We obtain below 
approximate solution using the WKBJ approximation. To implement the boundary condition at 
r = 0, under the WKBJ approximation, again it is better to transform to a new radial co-ordinate 
X, where r = e^. Substituting Hn{x) = /(a;) exp(— 3x/2)/KAr, equation (|50|) then becomes 




^3.V2 ( d\aNMn) ^ ^ (i(ajvM„) \ ^^^^ 
I dx^ dx I 

The boundary conditions now become f{x) — > as x ^ ±00. 

The WKBJ solution of the homogeneous equation for f = fn are 

/o r e2^\"^^^ /9 r e^^' \^^^ 
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Note that F^{x) satisfies the left boundary condition ^ as r = ^ 0, while 

F- staisfies the right boundary condition on /. The Wronskian of the two solutions is 

W = F^d{F_) / dx — d{F^)/dxF_ = —2. The particular solution to (^) , which satisfies the given 

boundary conditions can then be written as 



fix) 



dyG{x, y)^yl- { + 3 ^^";"^" ' 



where the Green function y) is given by 



1 



{B{y - x)F+ix)F^iy) + ^(x - y)F+(y)F_(x)] . 



(53) 



(54) 



Here 9{x) is the standard heavyside function, which is zero for x < and is unity for positive x. 
Transforming back to the original radial co-ordinate r we finally have 



Hpir,t) 



exp(r„t) 



dui-f^'^l [9{u - r)K{r;u) + 9{r - u)K{u;r)]x^^ { n^-^(a^M„(M)) 



1 d 



2KN{r) Jo r 
where we have defined 



V? du 



d 



du 



K{r;u) 



9 

4 + 2KN{r) 



-1/4 



4 2kn{u) 



-1/4 



exp 



'9 

4 + 2AtAr(r' 



r„r'2 V'^ dr' 



(55) 



(56) 



Of primary interest is the behaviour of Hp{0,t), the time evolution of the average current 
helicity associated with the fluctuating field, since it is this quantity which alters the alpha effect 
in the equation for the mean fleld evolution. This can be evaluated by putting r = in equation 
dH) . We get 



Hp{0,t) 



exp(r„t) 
V6 



dy- 



expi 



/_ ry di/_ [/£ I Kjv(0V^Ni/2 _ 3 
V Jo y' [U Kiv(y') ^ 2 



1 d 



(9 I «^iv(0)r \1 /4 

U + KNiy) I 



d aNM„{y). 



y^ dy V dy 2kn{0) 



(57) 



where we have used the identity (u/r)^/^ = exp(/"(3(ix/2x)) and also defined new variables 
y = u/a and y' = r' /a with = 2K7v(0)/r„. 

We can calculate Hp{0,t) explicitely once the functional forms of a]\f{y), AtAr(y) and Mn{y) 
are specified. For this let us consider the example of the previous section, where we TLL{r) is 
given by (^) , simulating Kolmogorov turbulence. In the previous section we saw that in the 
case Rm/Re » 1) all the modes were strongly peaked about a radius r = ^ lc{Rm/Re)~^^'^ ■ 

1 /2 

And the growth rate r„ ~ VI /I for a mode extending up to r ~ / with r„ ~ Re' V/L = for the 
fastest growing mode. As a model we therefore adopt a Gaussian form for M„(r) 



M„(r) = M„(0)exp[--^] = M„(0) exp(-6y2 



(58) 
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where b = {a'^/2rj) = {2r]/Tn){RrrJ RelD- Using the fact {2'q/Tc){Rm/ Rell) = 1, we have 
b = (Fc/r^) > 1 for ah modes. 



Further, since << Ic in general, in evaluating the integral in Eq. ( |5^ ) it will also suffice to 
take the form of T^Lif) in the region < r < / to get {K]\f{0)/Kj\f{y)) = (1 + 2by'^/3)~^. Also, we 
can model the scale dependent aAr(r) as follows: At any scale r > Ic, in the galactic context, the 
helicity structure function will scale as 

2C{0) - 2C{r) ~ {vf.Tr/h){nTr) = {n/hy, (59) 

where Vr and are the eddy velocity and correlation time at any scale r which staisfy the relation 
VrTr = r. Also C(r) is continuous at r = Ic, should have zero slope at the origin, and should vanish 
for r > L. A model helical correlation function which satisfies all the above requirements is given 

by 

2C(r) = ——[1 - for < r < L ; C(r) =0 for r > L. (60) 

So in the kinematic regime we have 

aiv(y) _ ^ 2 2 _ ^ 2 



2kn{0) 2KN{0)h'' Tnh 
Using the above forms of ctNiy), i^N{y) and Mn{y) in Eq. (57) we get 



y'a' = —y'. (61) 



Hp{Q,t) = Ci(6)-i^M„(0)exp(F„t) (62) 



where Ci is given by the integral 

ry dA/_ 



I foo ^ exp(_ Jq I4 -r i+2f)i/V.s^ 2 



,,'2 



C,{b) = ^l dy " ' ^ X 4^(/^(y^exp(-&y^))) (63) 

VeJo (| + _|i^)i/4 y^dyV dy'' 

The constant Ci can be comupted numerically for any b. We give below an approximate 
analytical estimate. First we note that the integrand is suppressed exponentially for y >> 1. It 
then suffices to expand the square root in the exponential, assuming y^/(l + 26y^/3) << 9/4. The 
exponential term then gives a term (1 + 26y^/3)~^/'^'' to the integrand of Eq. (63) . Consider first 



1 d f A d 



the case 6 = 1, appropriate for the fastest growing mode. After repeated integration by parts we 
then have 

5 eyi]i{-uP') 125 ex.p{—w'^'^ 



^-(1 + 10^2/9)5/4 ^"(l + 10u;V9)9/^ ^''^ 

Note firstly that Ci(l) > 0. Now y"e~^^ is strongly peaked about y = (n/2)"^/^ for large n One 
can then replace the remaining slowly varying parts in the integrands above by its value at the 
peak. The remaining integrals are gamma functions, which can be evaluated to give Ci ~ 0.05 for 
for the fastest growing mode. 
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1 /2 

For modes which extend up to r ~ L, with growth rate ~ F/L, 6 = Rj ~ 170. For 
b >> 1, 1/46 << 1 and the contribution to the integrand in ( p3| ) from the exponetial term is 
^{1 + 26yV3)-V4f' ^ 1. Also for h » 1, we have main contribution to the integral coming from 
regions ^ b^^; In this case an approximate evaluation of the integral along lines similar to (H) 



gives Ci(6) w 10-2^2, for large 6 >> 1. 

A number of comments are in order: Firstly, note that in the kinematic regime, H(0, t) grows 
exponentially. So the assumption of small-scale stationarity made for example by Gruzinov and 
Diamond (1994) to derive a constraint on H is not valid. Also since Ci is positive, the growth of 
small-scale fields and hence H{0,t) goes to decrease the effective alpha effect. The extent of the 
decrease depends on how much the small-scale field grows before it saturates. In the kinematic 
regime, a mode with growth rate r„ ~ vi/l, for example, leads to a reduction in the alpha effect 
to a value 



aeff = 2C{0) -4aH{0,t) 



h 



1 



4aCi{bn)Mn{0)e^"' 



r„L2 



2C(0) 



i-(4 



2 2Cl{bn) 



■ (65) 



Since Ci « 1 in general and Rarnbi{l) >> 1 from Eq. (48), we see that the reduction 
to the a-effect due to ambipolar drift is negligible even as the magnetic field energy grows 
to equipartition. A similar result obtains if we use the reduction in alpha effect estimated 
in the quasilinear approximation of Gruzinov and Diamond (1994). In our notation their 
a^ff = cecD = 2C(0) — TH{0,t)/{-Kpn), where r is a turbulence correlation time. Using H{0,t) 
determined above, one then gets for a mode with growth rate ~ vi/l, 



acD = 2C(0)[1 - i4Ci/3)fil,t){viTl)/L^ 



(66) 



Here we have assumed that the small scale field energy density is a fraction f{l,t) of the energy 
density associated with eddies of scale /. Again for Ci « 1, the correction to alpha in the 
quasilinear regime is modest, as long as f{l,t) < 1. Ofcourse, it is unclear what is the domain 
of validity of the quasilinear approach. One should take the estimate in Eq. (|66[), of alpha effect 
reduction at best as a rough guide, since neither we (nor Gruzinov and Diamond 1994) are not 
treating the full nonlinear MHD problem. Finally note that H = Ci{Vl/Tn)Mi/h « Mi/h, since 
(Q/Tn) « 1 and Ci << 1. So our neglect of the ajy coupling term for the evolution of M^, in 
section 5 is an even better approximation than was argued there. 



8. Non-linear effects due to current helicity 

As mentioned earlier the effect of ambipolar drift is to add to the diffusion coefficient an 
amount 2aML{0,t) and decrease by 4:aH{0,t), in both the mean field and the correlation 
function equations. We now consider some implications of the additional term due to the current 
helicity H{0,t). First let us look at the non linear evolution of H(r,t) itself. The general solution 
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of (^) for the evolution of H, when Ml is also settling towards a saturated state is beyond the 
scope of this work. The following general comments can however be made: the increase in kjv due 
to ambipolar diffusion will lead to a further damping of H. However the change (decrease) in a^v 
can lead in principle to a further increase in H due to the fact that Ml has a negative curvature 
near the origin. 

To illustrate this more interesting nonlinear effect, we consider the following simpler problem. 
We assume that Ml has attained a saturated state and also ignore the effects of turbulent diffusion 
for analytical tractability. The equation for H then becomes 

d - Iriamu-,^ (^'-pj]H = (r'-f{a{r)ML{r) - 4ai/(0, t)Mz.(r))) ^ p{r,t) (67) 

ot or \ or J r^ Or \ Or J 

where r]ambi = + 2aMi(0) as before and a{r) = 2C(0) — 2C(r). One can formally solve this 
equation using the Greens function for the operator on the LHS of Eq. (|67|). One gets 



rt /"OO /"CXD 

H{r,t)= dt' r'^dr'G{r,t;r',t')p{r',t') + r"^dr'G{r,t;r' ,0)H{r' ,0) (68) 
Jo Jo Jo 

The Greens function G can be obtained by standard methods (see for example Burton 1989); we 
get 

G{r,t;r ,t) =e{T)- T^^^M — ^ ] . (69) 

where we have defined t = t — t' and Is/2{^) = (2/(7ra;))^/^[coshx — sinhx/x] is the modified 
Bessel function. One can see that the Green function G is similar to that of the diffusion equation 
; infact it is akin to the radial Green function for diffusion in 5-spatial dimensions. We are again 
particularly interested in the evolution of the current helicity H{0, t). To look at its evolution take 



the limit r — > in Eq. (68) to get 



mt) = k-Y" dt' r'W0(T)— -^^exp[--A!^]p(r^^') 
3 vr Jo Jo {'iriambiTr' 8r]ambir 

+ k-)'" r'^dr'— exp[-/li^]i7(r',0) (70) 

3 vr 7o {,'^'nambitr'^ 8r]ambit- 

An initial distribution of H{r, 0) will be damped by diffusion. Also since the term proportional 
to a(r) in p(r, t) is constant in time the genetration of H{0, t) due to this term will also be 
damped eventually by diffusion. The only potentially interesting term for the further growth of 
H is the nonlinear term in p{r,t) oc H{0,t). To see its effect, let us assume that the steady state 
longitudinal magnetic correlation function can be described by ML(r) = Ml(0) exp[-rV(2i.^)]. 
We also take an intial H{r,0) = -ff(0, 0) exp[— r^/(2L|^)]. Using Eq. ( [70| ) , doing the r' integrals 
we then get 

_ 20aAfc(0) f _ g(0.n ^ 'iim pi) 

LI h (l+4%mWT/L?„)'/2 (l+4,„„nt/iH)''' 
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This is an integral equation for H(0,t). An approximate solution to this equation can be 
obtained as follows. Consider times < t < ti such that '^rjamhiti / < 1. For these times one 
may approximate the integrand in Eq. (^) by just H{0,t'). Then H(0,t) satisfies the simple 
differential equation 

' H{0,t) H{0,0) 



dH{0,t)/dt = 20aML(0) 



t2 t2 



(72) 



So if Lm < Lh, H{0, t) can grow exponentially. We saw in section 5 and 6 that Lm is generally the 
diffusive scale. Since during linear evolution of H{r,t), the current field correlations are generated 
from ML{r,t) through the alpha effect, we expect Lh ~ Lm- To estimate the maximum further 
growth of the current helicity, assume that we do have L^ < Lh, initially. Then for t < ti the 
current helicity grows exponentially with H{0,t) = if(0,0)[(l - fmn) exp{20aML{0)t/L^) + fmn], 
where fmH = Lm/Ln < 1- At time ti we then have H{0,ti) = H{0, 0)[{1 — fjnH)e^^'^ + fmn], where 
we have used the fact r]ambi = ?? + 2aML{0) ~ 2aML{0) . At later times, one can treat H{0, t') in 
the integrand of Eq. (^T|) as a slowly varying function compared to 1/(1 + 4:r]ambi{t — ^O/^m)"^^^) 
pull it out of the integral and integrate the resulting equation to get 



At large times one sees that H(0,t) — > 2H{0,ti) < 24ff(0,0). The upper limit obtains only 
when fmH « 1- A further growth of H{0,t) in the nonlinear stage implies a further reduction of 
cteff- Incase of purely ambipolar drift this reduction is still very small (cf. Eq. (|65|), while ugd 
reduction depends on how large f{l,t) grows before attaining the saturated state. For example if 
as argued in Paper I, f{l,t) « 1 in the saturated state then acD will also suffer only a modest 
reduction, due to small-scale dynamo action. 

Before ending this section, we mention two other interesting non linear effects, which obtains 
under the action of ambipolar drift, which both invlove the existence of a non-zero magnetic 
helicity: 



• (a) Relaxation through selective decay 

Suppose we start off with a random field configuration which has an initial non-zero magnetic 
helicity Pf/(0,ii), and ask how it will evolve under the action of ambipolar drift. Due to 
ambipolar drift and the resulting ion-neutral friction, the field energy is constantly drained 
into heat. This is refiected by the diffusion terms in the equations for Mi and H. However, 
for a large conductivity, the magnetic helicity is almost conserved, reflecting the fact that 
the field is almost frozen into the ions. A non-zero magnetic helicity also implies neccessarily 
a minimum non-zero field energy. So the random field cannot evolve into a zero energy 
configuration, but rather into a minimum energy configuration conserving magnetic helicity. 
Such a selective decay of energy, conserving helicity is thought to lead to self organisation 
into larger scale structures (Taylor 1974). Our formalism incorporating ambipolar drift offers 
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a simple route to study the dynamics of this relaxation process, and the final relaxed state. 
A study of such relaxation trough selective decay, using the equations for Mi and H derived 
here is in progress and will be reported elsewhere. 

• (b) Inverse cascade due to the non linear dynamo 

Suppose we had no kinetic helicity. But we had created large, random, small-scale fields, and 
also some small-scale magnetic (and current) helicity. For example, large random magnetic 
fields, correlated on small scales may be generated during phase transitions in the early 
universe. Further, in some of these phase transitions, like the electroweak one, there are 
speculations that large magnetic (and current) helicities may arise (Cornwall 1997, Joyce 
and Shaposhnikov 1997). Then our model non linear equation for the mean field (|l^), (and 
the work of Pouquet et al. 1976), indicates that the small-scale current helicity can lead to 
large-scale dynamo action. The alpha effect will be purely magnetic. This will lead to a 
coupling of small scale to large scale and " inverse cascade" (or dynamo growth) of magnetic 
energy to larger scales, due to purely nonlinear effects of the Lorentz force. If there were no 
other source of magnetic energy, the energy and current helicity would decay montonically. 
However the approximate conservation of magnetic helicty under near ideal MHD conditions, 
would still keep this dynamo active. It would be interesting to explore this issue further, 
using the full MHD equations, 

9. Discussion and conclusions 

We have revisited here the dynamics of fluctuating magnetic fields in turbulent fluids, ab 
initio. In doing this we have also incorporated the effects of ambipolar drift, as would obtain 
in a significantly neutral gas. Ambipolar drift introduces a magnetic field dependent addition 
to the velocity field in the induction equation. The resulting non linear equation may also be 
viewed, albiet with some caution, as a toy model for the MHD problem. Assuming that the 
velocity field has a turbulent component, we have derived the evolution equations for the mean 
and fluctuating magnetic field. These equations are used to discuss a number of astrophysically 
interesting problems. We summarise below the principle results of our work. 

First, in the presence of ambipolar drift, the dynamo equation for the mean field, and the 
equations for the magnetic correlations are both modified. Assuming a Gaussian closure, one gets 
an extra diffusion term proportional to the energy density in the fluctuating fields and a reduction 
to the alpha effect due to the average current helicity of the fluctuating fields. These equations 
( Eq. ( p!6[ ) , Eq. (^) and Eq. ( p5[ ) ) form a closed set of nonlinear partial differential equations 
for the evolution of both the mean magnetic field and the magnetic fluctuations, incorporating 
the back reaction effects of ambipolar drift. Due to the Gaussian closure approximation, the 
nonlinearity appears as time-dependent co-efficients involving only the average properties of the 
fluctuating field itself. 
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We applied our equations in section 5 and 6 to discuss small-scale dynamo action in galaxies, 

assuming a Kolmogorov type galactic turbulence. In the kinematic phase, dynamo action 
exponentiates small-scale fields provided the MRN associated with the turbulence, exceeds a 
critical value Rc ~ 60. Further, for Kolmogorov type turblulence, the critical MRN for excitation 
of a mode extending upto r ~ Zis also Rm{l) = Rc, as expected from the scale invariance in the 
inertial range. In galaxies, it is likely that the ratio of the magnetic to fluid reynolds number 
Rm/Re >> 1- In this case the fastest growing modes grow on a timcscalc comparable to the 

1/2 

turn-over time of the smallest eddies at the cut-off scale Ic, with T ~ VRj jL ^ Vc/lc- Modes 
whose longitudinal correlation function extends upto r ^ I > 1^ grow at a slower rate F ~ vi/l. 
Further, the field is strongly peaked about a region r = < lc{Rm/Re)~^^'^ about the origin 
for all the modes. For the most rapidly growing mode, w{r) changes sign accross r = Ic and 
rapidly decays with increasing r/lc- Interms of the Zeldovich rope-dynamo, one may picture, the 
small-scale field as being concentrated in ropy structures with thickness ~ and curved on a 
scale upto ~ Z for a mode extending to r ~ L 

As the small-scale fields grow, ambipolar drift adds to the diffusion coefficient and the effective 
MRN, Rambi for fluid motion on any scale of the turbulence decreases. If RambiiL) decreases to 
a value Rc ~ 60, dynamo action will stop completely. In a sufficiently weakly ionised medium 
ambipolar drift by itself can then lead to a saturation of small-scale fields before the turbulence 
is drastically modified. The required small ionisation fraction may obtain in the universe just 
after recombination and in proto stellar disks, but not in galactic gas. The effective magnetic 
Reynolds number in galactic gas, even including ambipolar diffusion, is much larger than Rc. 
In this case, as the the small-scale field grows in strength, it continues to be concentrated into 
thin ropy structures, as in the kinematic regime. In a companion paper (Subramanian, 1997 : 
Paper I) we have built upon the results obtained in sections 5 and 6, to discuss in detail how the 
small-scale dynamo may saturate in the galactic context, while preserving large-scale dynamo 
action. The crucial property of the small-scale dynamo generated field which allows this is its 
spatial intermittency. The field can build up locally to a level which will lead to small-scale 
dynamo saturation, while at the same time having a sub-equipartition average energy density, so 
that the diffusive property of the turbulence is not drastically affected. 

We also considered in sections 7 and 8 the evolution of the current-field correlations. The 
dynamics of helical correlations has not been as extensively studied in the literature as that of 
longitudinal correlations (and the small-scale dynamo) . In the kinematic phase the average current 
helicity associated with the small-scale field, H, decays due to diffusion, unless forced to grow by 
coupling to the exponentially growing longitudinal magnetic correlations. We showed in section 7, 
using a WKBJ approximation to the relevant Greens function, that the coupling with Mi^, leads to 
an exponential growth of the current helicity H{0, i) in the kinematic regime. So the assumption of 
small-scale stationarity made for example by Gruzinov and Diamond (1994) to derive a constraint 
on H is not valid. The growth of small-scale fields and hence -ff (0, t) goes to decrease the effective 
alpha effect. The extent of the decrease depends on how much the small-scale field grows before it 
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saturates. In general we find the reduction in the alpha effect due to the ambipolar drift itself, as 
given by Eq. (^) to be negligible in galaxies. The alpha effect reduction due to the growing small 
scale field, in the quasi linear approximation to the full MHD employed by Gruzinov and Diamond 
(1994) is also small (Eq. ( |66[ ) for aGo), as long as f{l,t) < 1. 

In section 8 we considered the non linear growth of the current helicity. We find that 
in principle H{0,t) can grow further, due to its nonlinear coupling with M^. This could 
lead to further reduction of acD- Even then, provided the small-scale dynamo saturates at 
sub-equipartition levels as argued in Paper I, the alpha effect is likely to be preserved until the 
large scale field itself grows to near equipartition level. Ofcourse, it is unclear what is the domain 
of validity of the quasilinear approach, and so this result must be treated with caution, until 
rederived in a full MHD treatment. 

The equations for the magnetic correlations incorporating ambipolar drift derived here, may 
offer a simple route to study the dynamics of relaxation through selective decay. In the case 
when magnetic helicity (without kinetic helicity) is initially present, non linear effects due to the 
Lorentz force can also lead to a magnetic alpha effect and dynamo generation of large scale fields, 
as envisaged in the early work of Pouquet et al. (1976). It will be worthwhile to study these issues 
further, possibly with a numerical treatment of the equations derived here. It is also desirable to 
return to the dynamics of small scale fields in the full MHD context; Perhaps by finding simple 
ways of incorporating the dynamics of the velocity correlations on an equal footing as that of the 
magnetic correlations. We hope to return to these issues in the future. 
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A. Evolution of fluctuating field correlations 

We give below more algebraic details involved in simplifying the first term in Eq. (|2l| ) . This 
term is given by 

^Rjpq (v^{y,t) ''Riim{vT{^,s)Mmq)^ ds >= -eituf-ulm.e^jrs(-spq Q^r Q^t '^^^ ^rnq (Al) 



For examining the evolution of Ml one needs to multiply the above equation by r'^r^ jr^ . We 
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can simplify the resulting equation by using the identity 

r r^Ti — TTi: = ^ ^ , - Sjtr'- dirV^ -—r - 6jtSirA (A2) 

where A = T^PMmq- Then using euueuim = ^iAm - ^imSu, and the definition of Tll,Tnn and C, 
straightforward algebra gives the contribution of the first term to {dMi/dt) 

^(1 St term) = -^U'tJ-^) + f + ,CH (A3) 
at or or 2 

The second term of Eq. ) gives an identical contribution. 



To derive the evolution of H due to these terms multiply Eq. ( [Al| ) by ^ijfr^ . Using the fact 
that the turbulent velocity and small scale field have vanishing divergence, we have Mj^ j = and 
Tjjj- = 0. This allows one to simplify the contribution from the first term to 

dH e-rrf 

— (1 St term) = [Tij,trMtr + TtrMij^tr - Tir,tMtj^r - Ttj^rMir^t] (A4) 

The first two terms on the RHS of Eq. ( |A4| ) can be further simplified by noting that eijfTij = 2Cr^ 
and EijfMij = 2Hrf . We have then 

- mj,trMtr + TM,,tr] = " [TllH" + T'^^H' + + MlC" + MiC' + ^^^] (A5) 

Here prime denotes a derivative with respect to r. To evaluate the contribution of the last two 
terms on the RHS of Eq. (|A4D it is convenient to split up the tensors Mjj and Tij into symmetric 
and antisymmetric parts (under the interchange of (ij) ). We put a superscript 5 on the symmetric 
part and A on the antisymmetric part. Then we can write after some algebra 



r 

Adding the contributions from Eq. ([A5|) and ([A6|) gives 



(A6) 



r)IT If) r) 

— (1 St term) = --j^ir'^mH + CM,] (A7) 



Qf > r'^ dr'' dr 

The second term of Eq. (^ ) gives an identical contribution. 



B. The WKBJ analysis of the kinematic small-scale dynamo 

First, in order to implement the boundary condition at r = 0, under WKBJ approximation, 
it is better to transform to a new radial co-ordinate x, where r = e^. Also to eliminate first 
derivative terms in the resulting equation we substitute ^{x) = exp{x/2)Q and get 

0+p(-)e = o (Bi) 
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where 



p{x) 



-{T + U)e 



2x 



1 

4 



The WKBJ solutions to this equation are hnear combinations of 

1 



1 

6 = exjp{±i J p^/'^dx) 



(B2) 



(B3) 



The solutions have to satisfy the boundary conditions Q{x) for x ±00. One therefore has a 
standard WKBJ eigenvalue problem for the determination of T. Note that asx— > — oo,p— > —9/4 
and so the WKBJ solutions are in the form of growing and decaying exponentials at this end. This 
also the case as x — > +00 since p — > — Fe^^ < for growing mode solutions with F > 0. In order 
to match the boundary conditions at both ends, the solution has to be the growing exponential 
at X = —00 and transit to the decaying exponential as x — > +00. This can only obtain if the 
p{x) goes through zeros, by U becoming negative for some range of r. For U considered here, in 
general, p{x) goes through two zeros. Suppose these occur at xi and X2 with xi < X2- Then the 
WKBJ soultions will be oscillatory in the range xi < x < X2- The requirement that the oscillatory 
solutions match on to the growing exponential near x = — cxd and the decaying exponential as 
X — > +00, gives the standard condition (cf. Jeffreys and Jeffreys 1966, Mestel & Subramanian 
1991) on the the eigenvalue F 

,1/2/ A _ (2n + l)7r 



X2 



p ' {x)dx 



(B4) 



One also determines the eigenfunction B(x), under the WKBJ approximation, to be 

A 



G 



(-p)l/4 

21/2A 



exp 



,1/4 



■ sm 



{-pf'^dx 
r{pfl^dx+^- 

Jx\ ^ 



(-1)"A 



exp 



{-pf^dx 



X2 



X < Xl 



Xl < X < X2 



X2 < X 



(B5) 



B.l. Critical MRN for the marginal mode 

Let us first use Eq. ( |B4| ) to determine the critical value of Rm = Rc needed for growth 
of the small-scale fields. For this we have to put F = and find Rm which satisfies ( p^ ). With 
F = 0, using U from ( ^ ) , defining y = r/L we have in the range Ic/ L < y < 1, 

_ -(9/4) - (29i?^/54)//3 + (13i?^/108)//3 

(l + y^/^Rm/Sy ^""^^ 

We will find that Rm = Rc is large enough that one can assume 1/Rm « V^^"^ in the above 
expression for p. In this case the zeros occur approximately at ri/L ~ (174/39i?c)^'''^ = 2/0 and at 
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r2/L ^ 1. (Here ri = exp(xi) and r2 = exp(x2) )• The integral condition then becomes 

174 



Jx, ^ ^ Jyo 12 39Rcy^/^^ y 2 



(B7) 



where we have looked for the critical reynolds number for the principal mode with n = 0. The 
integral in Eq. ( ) can be done exactly and gives the condition 



3 .13x1/2 
2^2 



ln\ 



1 + 



yl 



yo 



]-Vl-yo 



vr 
2 
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The solution of this equation implies a critical value for the excitation of the marginal mode 
Rm = Rc ^ 60. Note that Rc was estimated by taking n = in Eq. ( ). One can easily get 
also the limiting magnetic reynolds number needed for the exitation of higher order modes. It 
should also be pointed out that in Eq. ( ^ ) for p{y) , y^/^Rm = {r/L)'^/^Rm = {r/l)'^/^Rm{l) ■ 
So the above equations determening R^. are the same if we replace {L,Rm) by {l,R„i{l))- This 
shows that the critical MRN for excitation of a mode extendinng to r ~ / is also -Rm(0 = ^c, as 
expected from the scale invariance in the inertial range. 

The actual value of the MRN at the outer scale is likely to be much larger than Rc, in 
galaxies. Infact it is likely that Rm/Re >> 1- Let us now estimate the growth rate of the fastest 
growing mode in this case. 



B.2. Growth rate for the fastest growing mode 

For this we fix R^ and Re and look for the value of T which satisfies Eq. ( ) . For 
Rjn » -Re^the potential U is negative at r = L, decreases monotonically from r = L to r = Ic and 
is still negative at r = Ic- It only starts increasing for r < Ic- The fastest growing mode is then 
expected to concentrate at r < Ic- Also the turning point corresponding to x = xi occurs at r < 
for all the modes. So to determine the growth rate of the fastest growing mode and also examine 
the structure of the modes at small radius, one must adopt the form of T/^_L(r) with r < Ic, given 
in Eq. ( |3^ ). For r < Ic, one then has kjv = ^/[l + (l/3)(r/r(^)^] and 

V U-? \2z''^ - 1 - 4z2l 



Here z = r j (\/3rd) with 



Rm [Ic) 



setting the scale over which the potential varies. Using this form of the potential we then have 



Aqz^ - Bqz^ - 9/4 

P^'^ = (1 + ^2^2 (Bll) 
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~ 1/2 . — 

for the ranere 0<r <lc. Here T = r/{VRe' /3L) is a normalised growth rate, Aq = 15/4 — F and 
Bo = T- 1/2. 

For the above form of p{z), there is only one real positive zero zq. It turns out that zq is 
large enough such that the 9/4 in p{z) above can be neglected compared to the other terms 
giving 2:0 ~ (F — 1/2)^/^/(15/4 — f )^/^. Note that for the form of Tiiir) we have adopted 
the potential is discontinuous at r = Ic- For r = Ic — e, U = U~ = —AiVRj /3I/) while for 

1 /2 

r = lc + e,U = U+ = -A/'i{VRe' I'iL). As we mentioned earlier this does not alter the results 
qualitatively since they are generally dependent on integrals over U . It means however that for 
4/3 < F < 4, which we will see obtains for the fastest growing mode, the outer zero of p is at 
r = Ic OT z = Zc = Ic/iV^Td) = (i?m/3i2e)^''^- For f corresponding to the fastest growing mode, 
we then obtain the integral condition 

[Aoz^-Bo^dz ^ vr ^^^^^ 

Zq z'^ Z 1 

In the above we have again assumed that Zq » 1, which we will show to be true below. This 
integral can be done exactly and gives the condition 

Al'hn[2AQzl/BQ - 1] - 2{AqzI - Bq)/zc = |. (B13) 

In the above condition, since z^ = {R^/'iRe) » !> one can get a good iterative approximate 
solution for f . In the first approximation one neglects the constant terms compared to the z^ 
terms to get 

f ^ 15/4 - {^/2ln{Rm/Re)f. (B14) 

For example for Rm/Re = 10^^ one has f ~ 3.748 or F ~ So as advertised the 

fastest growing modes grow on a timescale comparable to the turn-over time of the smalles, cut-off 

1/2 

scale eddies, if Rm » Re, with F ~ VRe /L. One can also go back and check, using this value 
of F that ~ 30 >> 1 and so the approximations made assuming zo,Zc » 1 are very good. 



B.3. Spatial structure of the eigen modes 

We now briefly consider the spatial structure for various eigenmodes of the small-scale 
dynamo, for the case Rm{lc) = Rm/Re » 1- In this case as we mentioned earlier, the turning 
point corresponding to x = xi occurs at r < Ic for all the modes. The eigenfunction for r < Ic are 
then given by Eq. ( |B5| ) with p given by the form in ( Bll| ) . 



Infact near the origin, one can use the the original equation (^) to find the behavior of Ml 
and w{r). One finds that for small r, <I> satisfies the equation 

2^> 

^-^ + o4. = (B15) 
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where the constant a = (5 — F). The solution can be found by elementary methods 
$ = z'^sm{^/az)/{^/az). So then have 



ML{r,t) 



ML{0,t) sin(^2:) 

[1+^2] 1/2 



(B16) 



A similar result can also be found from the WKBJ solution, where one gets for z < zq, 

e^*F{z) 



with 



F{z) 



Ml 



A, 



j3/2(l + z2)l/2 



{-P{Z)Y" 



exp 



(-p(z))V2^' 

20 Z 



As z ^ 0, —p{z) ^ 9/4 — az'^ and so this soln goes over to 



Ml 



ML{0,t) 



1 



[l + z2]l/2 (1- (4q/9)z2; 



■exp[ 



-az 



(B17) 

(B18) 

(B19) 



This is in good agreement with the exact soln determined for small z. For larger z, Ml decreases 
monotonically with z. The WKBJ treatment gives for zq < z < Z2 {z2 is the outer turning point) 



Fiz) 



(p(z))l/4 



sm 



(B20) 



Away from z = zq, the z'^ term in p{z) dominates and p (15/4 — F) = Aq. So for z » zq we 
have 



Ml ~ — rn sm 



(B21) 



The above equations show that Ml decreases rapidly with increasing z. Since z = r/r^, 
one sees that Ml and hence w{r) for all the modes are strongly peaked about the radius 

1/2 

r ^ rd = Ic/Rm {Ic) « Ic, for the case Rm{lc) » 1- 

We also mentioned that w{r) must become negative at some radius. Prom the WKBJ solution 
it is apparent that the number of zero crossings for the WKB solution will depend on the order 
n of the mode. Let us consider the fastest growing mode. We saw earlier that for this mode the 
outer turning point is at r = /c or z = Zc. For this mode one can check that from the WKBJ 
solution that w{r) > 0, in the region < r < Ic- For r > Ic and Rm/Re >> li one has from Eq. ( 
-{AVRl^'^/9L){r/lcy^/^ and kn ~ {VL/3Re){r/lc)^/^. Therefore 

F 



m ) . U{r) 



for lc<r<L; R^/R, » 1 



(B22) 



For the fastest growing mode one can igmore the constant term and get —p{r) ^ T{r /Ic)'^^^ . Then 



ML{r) oc i-y/'^exp -3f^/\{r/k)^/^ - 1) 



(B23) 
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and 

1 f2,lc. 



w tx 



(^)l/3_pl/2 _3f l/2((^//jl/3 _ (B24) 



y2 V3 V 

For r ~ 3.75 we can see that w^r) < for r > 1^ and its modulus decreases to zero rapidly with 
increasing r/lc- So w changes sign accross the transition point r = Ic for the form of longitudinal 
correlation function T^l we have adopted. For modes with a smaller growth rate F, we see from 
Eq. ( B22| ), that the outer turning point, got by putting j5(r) = 0, is at r/l ^ (13/(12?;))^/^, where 



F; = T /[vi/2>l). So modes with growth rate F ~ vi/l, extend upto r ^ I. 

In summary, in the case Rm/Re » 1) w{r) is strongly peaked about a region 
r < lc{Rm/Re)~^^'^ about the origin for all the modes. For the most rapidly growing 
mode, w{r) changes sign accross r = Ic and rapidly decays with increasing r/lc- (If one had 
adopted a sufficiently smooth form for Tn^r) around r = Ic, the value where w changes sign would 
still be r ~ Ic, but could have been better determined by the WKBJ analysis). Also slower growing 
modes with F ~ vi/l, extend upto r ~ L A more thorough analysis of the eigenfunctions can be 
found in Kleeorin et. al. (1986), for the simple case when the longitudinal velocity correlation 
function has only a single scale. Their analysis is also applicable to the mode near the cut-off scale 
of Kolmogorov type turbulence. 

Let us now consider the corresponding eigenfunction for the marginal mode. In this case, as 
we saw earlier, the inner turning point occurs at r = ri = Li/q ~ 0.14L and the outer turning 
point is at r = r2 ~ L. The WKBJ solution can be used to determine the the eigenfunction. We 
get for r < ri, 

where p{y) is given by Eq. ( ^6| ) . For r « L, we can neglect the y^^^ term in have Eq. ( ^6| 
) compared to the constant and y'^/^ terms. This implies Q{y) ~ y^^^ exp(— 13i?my^'''^/54) and so 
for the marginal mode 

M ( f\ exp(-(13/54)i;^y^/3) 

(1 + (l/3)i?„,V3)i/2 (B26) 

One sees that the eigenfunction is concentrated in a radius r ~ L/Rc = rc for the marginal 
mode. Infact if we define z'^ = Rcy^^^ /3, then for the marginal mode with F = we have from Eq. 
( ), p = (-9/4 - 25zVl8 + 13z'^/12)/(l + z'^)^. This is very similar to the p{z) defined in Eq. 



Bll| ) except for the identifications z ^ z, Aq ^ 13/12, Bq 25/18. Going through the same 



analysis as for the p{z) of Eq. ( [Bll| ) , gives the properties of of w{r) for the marginal mode : For 
the marginal mode w{r) peaks within z ~ 1, corresponding to a radius r ~ L/Rc^^, changes sign 
to become negative at r ~ L and dies rapidly for larger r/L. 
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